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The explicit expressions for the non-singlet DIS structure function gi at small x are obtained by 
rcsumming the leading logarithmic contributions. The role played by the fits for the initial parton 
densities currently used in the DGLAP on the small- a; behavior of the non-singlet gi is discussed. 
Explicit expressions combining DGLAP with our results are presented. 
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The non-singlet component of the spin structure function g\ have been investigated in great detail in deep inelastic 
H 1 scattering (DIS) experiments. The standard theoretical framework for studying the DIS structure functions is provided 
by DGLAPp]]. In this approach, g^ s {x,Q 2 ) can be represented as a convolution of the coefficient functions and the 
evolved quark distributions. Combining these results with appropriate fits for the initial quark distributions, provides 
a good agreement with the available experimental data. 

However, the DGLAP evolution eqs. were originally applied in a range of large x values, where higher-loop contri- 
butions to the coefficient functions and the anomalous dimensions are small. Such corrections are becoming essential 
ly-j \ when x is decreasing, so DGLAP should not work so well at x <§; 1. Nevertheless, DGLAP predictions are in a good 
agreement with available experimental data. It leads to the conclusion that the impact of the higher-order corrections 
is negligibly small for the available values of x. Below we use our results |2( to show that the impact of the high-order 
corrections on the Q 2 and x -evolutions of the non-singlet structure functions is quite sizable and bounds the region 
of strict applicability of DGLAP to x > 10~ 2 . We also show that the reason for the success of DGLAP at x < 10~ 2 
Oh, is related to the sharp x -dependence assumed for the initial parton densities, which is able to mimic the role of 
high-order corrections. 

The paper is organized as follows: In Sect. 2 we discuss the difference of our approach with DGLAP. Then we 
compare our and the DGLAP formulae for asymptotics of g\. In Sect. 3 we suggest a method to combine DGLAP 
with our approach in order to obtain equally correct expressions for both large and small values of x. Sect. 4 contains 
our conclusions. 



II. COMPARISON OF DGLAP AND OUR APPROACH 

As the DGLAP -expressions for the non-singlet structure functions are well-known, we discuss them briefly only. 
In this approach, 9i%clap( x > Q 2 ) can be represented as a convolution 

9idglap(^ Q 2 ) = / -C(x/y) Aq(y, Q 2 ) (1) 

j x v 

of the coefficient functions C(x) and the evolved quark distributions Aq(x, Q 2 )- Similarly, Aq(x, Q 2 ) can be expressed 
through the convolution of the splitting functions and the initial quark densities 8q(x w 1, Q 2 ~ M 2 ) where /i 2 is 
the starting point of the Q 2 -evolution. It is convenient to represent f(x,Q 2 ) in the integral form, using the Mellin 
transform: 



diii r dh 2 i 

9i dglap( x i Q 2 ) = «/2) / — (l/xrC(w)6q(u;)e XP 7 ( W ) / -f a s (k 2 ± ) (2) 
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where C(u>) are the non-singlet coefficient functions, 7(0;) the non-singlet anomalous dimensions and Sq(uj) the Mellin 
transforms of the initial non-singlet quark densities. The standard DGLAP fits Sq(x) for the non-singlet parton 
densities (see e.g. Refs. 00) consist of the terms singular when x — > and the regular in x part. For example, the 
fit A of Rcf. 3] is chosen as follows: 



Sq(x) = Nrix~ a c/)(x) , (3) 
<t>{x) = (l-xf{l + fX S ) , 

with N, 77 being the normalization, a = 0.576, (3 = 2.67, 7 = 34.36 and 5 — 0.75. As the term x~ a in the rhs of 
Eq. is singular when x — > whereas the second one, <f>(x) is regular, we will address them as the singular and 
regular parts of the fit respectively. Obviously, in the to -space Eq. (|2J| is a sum of the pole contributions: 

00 

Sq(w) = Nt} (lo - a)^ 1 + ^2 m k f(w + k - a) -1 + 7(cj + k + 1 - a) -1 ) , (4) 

k=l 

with nik = l)..(/3 — k+ l)/kl, so that the first term in Eq. (@J (the leading pole) corresponds to the singular term 

x~ a of Eq. © and the second term, i.e. the sum of the poles, corresponds to the interference between the singular and 
regular terms. In contrast to the leading pole position to = a, all other poles in Eq. J3J have negative values because 
k — a > 0. An alternative approach was used in Refs. [? ], by introducing and solving infrared evolution equations 
with fixed a s . This approach was improved in Refs. Q, where single- logarithmic contributions were also accounted 
for and the QCD coupling was regarded as running in all Feynman graphs contributing to the non-singlet structure 
functions. In contrast to the DGLAP parametrization a s = a s (k 2 _), we used in Refs. another parametrization 
where the argument of a s in the quark ladders is given by the time-like virtualities of the intermediate gluons. Refs. 
suggest the following formulae for the non-singlet structure functions: 

/"IOC J 

g? s (x, Q 2 ) = (e 2 g /2) / £L(l/xrC NS (u>)8q(u) exp (H NS (u;)y) , (5) 
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with y — ln(Q 2 //i 2 ) so that /j, 2 is the starting point of the Q 2 -evolution. The new coefficient function Cats are 
expressed in terms of new anomalous dimensions H^s whereas if/vs account for the total resummation of the double- 
and single- logarithmic contributions (see Ref. 0] for details). 



III. COMPARISON OF DGLAP AND OUR SMALL-i ASYMPTOTICS 

When x — > 0, one can use the saddle point method in order to estimate the integrals in Eq. (0 and derive much 
simpler expressions for the non-singlet structure functions: 

5 f s -e^Mr , (6) 

with £ = Q 2 / (x 2 fi 2 ) and with the intercept loq = 0.42. Eq. JJjJ predicts the asymptotic scaling for the non-singlet 
structure functions: Asymptotically, s depends on one argument £ instead of depending on x and Q 2 separately. 

When the standard DGLAP fits, e.g. the fit of Eq. ©, are used, the asymptotics of g^DCLAP^^Q 2 ) ^ s a ^ so ^ ne 
Regge-like: 

(\ 7(a) lb 
(ln(Q 2 /A 2 ))/(ln( M 2 /A 2 ))) , (7) 

with 6= (33- 2n/)/127r. 

Comparison of Eq. ©and Eq. (JJJ demonstrates that both DGLAP and our approach lead to the Regge asymptotic 
behavior in x. However, it is important that our intercept loq is obtained by the total resummation of the leading 
logarithmic contributions and without any assumption about fits for Sq whereas the DGLAP intercept a in Eq. is 
generated by the phenomenological factor x~ ' 57 of Eq. J2J which mimics the total resummation. In other words, the 
impact of the higher-loop radiative corrections on the small- ir behavior of the non-singlets is, actually, incorporated 
into DGLAP phenomenologically, through the fits. It means that the singular factors can be dropped from such fits 
when the coefficient function includes the total resummation of the leading logarithms and therefore in this case fits 
for Sq can be chosen as regular functions of x. 
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IV. COMBINING DGLAP WITH OUR HIGHER-LOOP CONTRIBUTIONS 



Eq. accounts for the resummation of the double- and single logarithmic contributions to the non-singlet anoma- 
lous dimensions and the coefficient functions that are leading when x is small. However, the method we have used does 
not allow us to account for other contributions which can be neglected for x small but become quite important when 
x is not far from 1. On the other hand, such contributions are naturally included in DGLAP, where the non-singlet 
coefficient function C dglap and anomalous dimension "f dglap are known with the two-loop accuracy: 



C DGLAP = l + ^^C^\ (8) 

idglap - -^^r + {-^^) r 



Therefore, we can borrow from the DGLAP formulae the contributions which are missing in Eq. (jSJ by adding 
C dglap and 7 dglap to the coefficient function and anomalous dimension of Eq. (JSJ. ft is important to avoid a 
double counting DL and SL terms common for these expressions. 

In order to do so, let us consider the region of x ~ 1 where the effective values of w in Eqs. ll'lot are large. In this 
region we can expand -ff/vs and Cms into a series in l/u>. Retaining the first two terms in each series, we arrive at 
C NS = C NS + 0(a 2 s ), H NS = H NS + 0(a 3 s ), with (see Ref. for details) 



Cw^ + ^fiy + i/a;], (9) 

Hns = [2/0, + I] + (^E)\l/ U ) [2/- + I] 2 + D[l/ U + 1/2]. 

Now let us define the new coefficient functions Cms and new anomalous dimensions Cms as follows: 
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These new, "synthetic" coefficient functions and anomalous dimensions of Eq. (|l()|l include both the total resumma- 
tion of the leading contributions and the DGLAP expressions in which a s (Q 2 ) is replaced by A(ui) defined in Refs. 
because the factorization of the phase space into transverse and longitudinal spaces used in DGLAP to parametrize 
alpha s is a good approximation for large x only. 



[1] G. Altarelli and G. Parisi, Nucl. Phys., B126 (1977) 297; V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. textbfl5 

(1972) 438; L.N. Lipatov, Sov. J. Nucl. Phys. 20 (1972) 95; Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 
[2] B.I. Ermolaev, M. Greco and S.I. Troyan, Nucl.Phys. B594B (2001)71; ibid B571(2000)137; Phys. Lett. B579, 321,(2004); 



hep-ph/0503019 
G. Altarelli, R.E 


. Ball, S. 


A. Vogt. hep-ph/0408244 



[3] G. Altarelli, R.D. Ball, S. Forte and G. Ridolfi, Nucl. Phys. B496 (1997) 337; Acta Phys. Polon. B29(1998)1145; 

[4] 



